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We consider quantum theoretical effects of the sudden change of the bound- 
ary conditions which mimics the occurrence of naked singularities. For a sim- 
ple demonstration, we study a massless scalar field in (1 + l)-dimensional 
Minkowski spacetime with finite spatial interval. We calculate the vacuum 
expectation value of the energy-momentum tensor and explicitly show that 
^ ■ singular wave or thunderbolt appears along the Cauchy horizon. The thunder- 

bolt possibly destroys the Cauchy horizon if its backreaction on the geometry 
is taken into account, leading to quantum restoration of the global hyperbolic- 
ity. The result of the present work may also apply to the situation that a closed 
string freely oscillating is traveling to a brane and changes itself to an open 
string pinned-down by the ends satisfying the Dirichlet boundary conditions 
on the brane. 
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1 Introduction 



A question of whether the cosmic censorship holds is one of the most important issues in 
classical general relativity, but still remains far from being settled. Many attempts have so 
far been made to prove the validity of this conjecture and have achieved partial success under 
some particular setup. On the other hand, some of those attempts have resulted in rather 
showing the counter examples, indicating that the classical framework of general relativity 
admits appearance of naked singularities (See e.g., Review [1]). In the present situation of 
research on the cosmic censorship hypothesis described above, it is probably more sensible to 
classify naked singularities and ask whether it is a disastrous one for physics or a practically 
harmless one. 

In the previous work [2], as a physically sensible classification of naked singularities, we pro- 
posed the wave approach to probe naked singularities. There, we studied the wave propagation 
in various static spacetimes with naked singularities. The essential point was the uniqueness 
of self-adjoint extension of the time-translation operator A, defined by the spatial derivatives 
in the wave equation. We call the classical singularity wave-regular if the self-adjoint exten- 
sion is unique, i.e., there is no ambiguity in the boundary condition at the singularity. In the 
wave-regular case, thus there exists a unique solution to the wave equation for a given initial 
value; the predictability holds at least for the probe fields even in the classical framework. 
In this sense, the concept of wave-regularity is closely related to Clarke's idea of generalised 
hyperbolicity [3]. 3 Such a would-be singularity can be regarded as a harmless singularity. On 
the other hand, if the time evolution of an initial value is not unique so that the predictability 
for a probe wave field breaks down, such a naked singularity is called wave- singular. 

The idea of the wave probe rather than particle probe is based on semi-classical theory, in 
which geometry is treated classically but matter fields quantum mechanically, with the probe 
waves identified with the mode functions for quantum fields. One might anticipate that, once 
a theory of quantum gravity is established, one would not need be afraid of the emergence of 
naked singularities; a naked singularity could be resolved by replacing it with something like a 
boundary of Planck energy scale, and in the neighborhood of the resolved singular boundary, 
the predictability of physics would be restored quantum mechanically. But even in such a case, 
a remaining question will be how such an extremely high energy boundary looks like from 
distant observers who live in a low energy scale governed by (semi-) classical theories. Then, it 
seems plausible that a naked 'singularity' resolved by any kind of quantum gravitational theory 
looks wave-regular from distant observers, as the wave-regularity respects the predictability of 
low energy (semi-) classical physics. In this sense, naked singularities can be classified into 
two types: one of which can be resolved to be wave-regular by quantum gravity, while the 
other should be somehow prohibited. 4 Here, we would like to propose a version of (quantum-) 
cosmic censorship advocating that, in semi-classical theory, wave-singular naked singularities 
are prohibited to appear by quantum field theoretical effects. 

In the present work, aiming at examining the above version of quantum cosmic censorship, 
we shall focus our attention on the wave-singular case in which a timelike singularity emerges 

3 The question of hyperbolicity has also been discussed in conical spacetimes [4], and in spacetimes with 
hypersurface singularities [5] . 

4 Recently, this kind of classification of singularities has also been argued from the stringy theoretical view 
point (See e.g., Review [6]). 
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at some stage and the boundary condition is not a priori specified. This feature may manifest 
itself when the singularity appears in occasion of gravitational collapse. A simplified version will 
be the case that the boundary condition suddenly changes at some time so that the boundary 
becomes a naked timelike singularity. 

We shall consider such a case in the context of quantum field theory. The time-translation 
operator A for a field changes in time by the change of its domain T>(A) so that we naturally 
expect particle creations from a vacuum state of the field. We shall first show that this is indeed 
the case in the perhaps simplest example: a massless scalar field theory in a (1 + 1) -dimensional 
spacetime R x [0, 7r], with the boundary condition at and n suddenly changing from the 
Neumann to the Dirichlet boundary condition. We shall compute the vacuum expectation 
value of the energy-momentum tensor and see that it just gives the standard Casimir energy 
inside the domain of dependence while it is divergent at the Cauchy horizon for the initial 
hypersurface {0} x [0,7r]. We will observe the thunderbolt effect [7] in the sense that there 
appears a delta function like singularity which has a support only on the Cauchy horizon. If 
the backreaction on the geometry is taken into account, the thunderbolt effects will destroy the 
Cauchy horizon, resulting in quantum mechanical restoration of the global hyperbolicity. We 
shall also give discussion for the general boundary conditions case in the appendix. 

Long ago Anderson and DeWitt [8] discussed a similar problem in the context of topology 
change in spacetime which exhibits a splitting of the compact universe like 'trousers.' They 
discussed particle creations with an infinite amount of energy which appear an infinitely bright 
flash emanating from the crotch of the trousers. In a sense the emergence of a naked singularity 
is a kind of topology change which alters the domain of definition for quantum fields. Not 
surprisingly their result, which is not fully quantitatively described, is qualitatively very similar 
to ours. 

It is interesting to point out a possible application to string and brane theories [9]. In 
the brane picture of universe, an open string can propagate along a brane with both the ends 
satisfying the Dirichlet boundary conditions on the brane but also can pinch off from the brane 
to the bulk spacetime as a closed string. In that process the boundary condition at the ends 
changes. The emerged closed string will be highly excited because of the sudden change of the 
boundary condition as an analog of the phenomenon described in the present work. 



2 Hair of wave- singular naked singularity 

For self-containedness we briefly recapitulate the self-adjoint extension of the time-translation 
operator in a static spacetime case. The relevant mathematical materials are collected in the 
appendix of Ref. [2]. For thorough study of this issue, see Ref. [10]. 

When probing timelike singularities in a static spacetime with a scalar wave 0, the wave 
equation can take the form 

d 2 

^0 = -40, (1) 

where t is the Killing time and A denotes an operator containing only spatial coordinates and 
derivatives. We will find the preliminary time-translation operator A as a symmetric operator 
on a given Hilbert space 7i with domain C^°, a set of smooth functions with compact support 
on t = const, surface. A trivial symmetric extension is immediately made by taking its closure. 
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Consider then the sets IC± := {0± G H\A*(p± = ±i(p±}, which are called the deficiency subspaces 
of A, and the pair of numbers := (dim/C + , dim/C_) called the deficiency indices of A. 

If A with domain V(A) is a closed symmetric operator with n + = n_, then A has self-adjoint 
extensions. Let U be the partial isometries /C+ — > /C_. Then the self- adjoint extensions Ae can 
be obtained by taking the domain as [10] 

V(A E ) : = {0 O + + + [70+ | O e 0+ G /C + }. (2) 

Given self-adjoint extension Ae, dynamics of the scalar field for any initial data (0(0), 9 t 0(O)) G 
H x H can be defined as 




The point is that, once a self-adjoint extension is given, the vector 0(t) G 7i is defined not only 
inside the future domain of dependence for the initial hypersurface but also even outside the 
domain of dependence [11]. 

In general, self-adjoint extension is not necessarily unique, when spacetime fails to be glob- 
ally hyperbolic. Even when spacetime admits naked singularities, if n + = n_ = 0, then A has a 
unique self-adjoint extension and dynamics of the probe field is uniquely determined everywhere 
from the initial data by Eq. (3). This case is referred to as wave-regular. If n + = n_ = N ^ 0, 
then the partial isometry U is represented by an N x N unitary matrix U(N), and A has 
infinitely many different self-adjoint extensions, which are in one-to-one correspondence with 
U(N). In this case the naked singularity is called wave-singular. Each self-adjoint extension 
corresponds to the different boundary condition at the singularity and accordingly describes 
different time evolution of the wave for the same initial data. In other words, a wave-singular 
naked singularity has degrees of freedom for the possible choice of the time-translation operator 
A E . The degrees of freedom can be interpreted as the character or the hair of the wave-singular 
naked singularity, described by the isometries U(N). For example, when a scalar probe field is 
considered in (l + l)-dimensions of the finite spatial interval with the wave-singular singularities 
at both the ends, the hair of the singularities is U(2). 

It should be commented that the wave-regularity depends on the choice of our Hilbert space 
H. As emphasized in Ref. [2], the Sobolev space, which requires L 2 of the derivatives of as 
well as itself, will provide a physically sensible Hilbert space for the singularity probe. In 
the present work, however, as a simple illustration, we shall take L 2 -space of square integrable 
functions as adopted in Ref. [12]. The analysis in the next section can be straightforwardly 
generalized to the Sobolev space case. The results will somewhat change depending on the 
choice of the Hilbert space. 

3 Quantum behavior of a probe field 
3.1 A simple model 

In this section, as one of the simplest cases, we shall consider a quantum field theory of a massless 
scalar field in (1 + l)-dimensional Minkowski spacetime R x [0, it] 3 (t, x). As mentioned before, 
the hair for the boundary "singularities" in this case becomes the group U(2), but here we shall 
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consider only its subgroup U (1) x U(l) with each U (1) corresponding to the arbitrariness of the 
boundary conditions at x — 0, n, respectively. We are concerned with quantum field theoretical 
effects in the situation that the hair U(l) x U(l) suddenly changes at some time, say t — 0, 
which can be interpreted as a sudden appearance of wave-singular timelike singularities at 
x = 0, 7r for t > 0. 

The scalar field obeys the Klein-Gordon equation (1) with an elliptic operator 

d 2 

(4) 



dx 2 



We require the square integrability for the scalar field, i.e., G L 2 (0, 7r) = {0 | J* dx\4>\ 2 < oo }, 
so that we can regard A with V(A) = C£°(0, n) as a positive symmetric operator in the Hilbert 
space L 2 (0,7r) with the inner product, 

2 ) : = / dx 0*(x)0 2 (:r) . (5) 



Then, having a suitable self-adjoint extension Ae of A, we can obtain the solution for any initial 
data (0(0, x), 9*0(0, x)) from Eq. (3). 

We shall focus on a particular situation that the boundary conditions at x = 0, 7r are the 
Neumann boundary conditions, d<p/dx = 0, before t — but switch to the Dirichlet boundary 
condition, = at a; = 0, n, after £ = 0. This sudden change of the boundary conditions can 
be viewed as a particular case of the occurrence of wave-singular naked singularities at x = 0, re 
after t — 0. 



3.2 Mode functions 

Before t — the positive frequency mode functions satisfying the Neumann boundary condition 
at x = 0, 7r are given by 



fi m \t,x) 



-int 



cosnx , (n — 1, 2, • • •) , 



(6) 



which are normalized by the Klein-Gordon inner product: % {(/, c^g) — (9</, g)}. Then we can 
expand the field (f>(t, x) for t < as 



oo 

0(t,x) = /„(*,*) + ^ {af)/(-)(i,x) + a^fi in) *(t,x)} , 



(7) 



n=l 



where 



fo(t,x) = —={q + pt) , 
V 71 " 



(8) 



is the zero-mode. The expansion coefficients are required to satisfy the canonical commutation 
relations 



Jin) (m)t 



(in) (in) 
u n > a n' 



„(in)t (in)t 
a n 5 "n' 



, [q, p] 



(9) 
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and a£ n \ d£ n ^ commute with q and p so that the canonical quantization is reproduced: 

[<f)(t, x), d t (f)(t, x 1 )} = iS(x - x 1 ) . (10) 

Now let us consider the time evolution of <f>(t, x) into the future (t > 0) region of the 
initial surface t — 0, where the probe field (p{t, x) is supposed to satisfy the Dirichlet boundary- 
condition (p(t, 0) = (p(t, 7r) = 0, and accordingly the positive frequency mode functions are given 
by 

„— imt 

ft t) (t,x) = ^smmx. (11) 



In order to study the time evolution of (p into t > region, we need to modify the initial data 
at t = so as to fit the Dirichlet boundary condition. We shall discuss the positive frequency 
modes fn^ and the zero-mode fo separately. 
First, using the Fourier expansion, 

"2 ^ / 1 I 



cos //.(■= ./ — } ( 1 ) sin ma;, (12) 

ii \m + n m — n I 

m=l 



(m±n:odd) 

we can express the initial data of the positive frequency modes as 



f!T>(0,x) = f; (^— + —) smmx, (13) 

71 V n ^-^ \m + n m — nj 

m— 1 x / 

(m±n:odd) 

^fi in \0,x) = -infW(0,x). (14) 

The evolution of the mode function into t > region is obtained by inserting the above 
initial data (13) and (14) into Eq. (3), in which V(A E ) is now spanned by the basis {/i™ 1 '} so 
as to respect the Dirichlet boundary condition. Then we can find 

oo 

ft\t,x) = J2{ a ^ft t] (t,x)+(3 nm ft t] \t,x)} , (15) 

m=l 

with the coefficients 

2 pm 1 _2 frn 1 

C^nm \ Om±n:odd ; Pnm \ , Qm±n:odd i \ ®) 

7T V n m — n tt \ n m + n 

where and hereafter 5 m ± n:oc id means 5 m ± n:oc id = 1 if m ± n = an odd number and S m ± n]0 dd = 
otherwise. 

The consistency of the canonical commutation relations (9) requires the unitarity relation 

oo 

^2 (®nm®* nm , ~ PnmPnm') = tm' • (17) 

71=1 
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This can be checked by using the formula 



£ 

n = l 

(n:odd) 



n 2 ~ 8 



(18) 



Next, let us expand the mode which is the time evolution after t — of the zero- mode 
before t = in terms of the late time mode functions (11) as 



m=l 



Mix) = {c m fkr t \t,x)+cyw*(t,x)} 

The initial data is then given by 

m=l v v 



(19) 



d_ 
dt 



m=l V V 



(20) 
(21) 



by continuity at t — 0. Then, by using the formula: (sin mx, sin nx) = |5 m>n , the expansion 
coefficients are obtained as 



2 1/ i 
C m = — ^= g H p ] <5 m:0( ta . 



(22) 



3.3 Particle creations 

Suppose the quantum field (f> is expanded in terms of the (out)-mode functions as 
<f>(t,x) = ^{a^f^\t,x) + ^f^\t,x)} . 

m 

Then, from Eqs. (7), (15), and (19), we have the Bogoliubov transformation 

(out) n , V^f M_i_/5 flO) 1 "! 

~ / ^ | u -nm u ' n ~ ynm^n ( ; 
n 



(23) 



(24) 
(25) 



This clearly exhibits the mode mixing of the positive and the negative frequency parts, which 
implies particle creations in quantum mechanics [13]. 

The expansion coefficients of the (out)-mode functions are now required to satisfy the canon- 
ical commutation relations 



(out) (oirt)t 



(out) (out) 
u m 5 u rn' 



' (out)\ (out)t 



o, 



(26) 



7 



so that am and a™ are an annihilation and a creation operator, respectively, as well as 

and an™"* . Note that, for the unitarity, the contribution from the zero- mode is indispensable. 
Then, defining the (in)-vacuum state as 

af ) |0 i „>=p|0 m >=0, (27) 

we can construct the Hilbert space of quantum states in the Fock representation. 

Now we can evaluate the expectation value of the number operator N m : = a&T^ at"*' of 
the m-mode particle in | in > as 

n 

4 f <0 in \q 2 \0 in > 1\ 1 

This is finite, but the total number of the created particles J2m=i < ®in\N m \ 0i n > is logarith- 
mically divergent, so that the Bogoliubov transformation is not unitarily implementable [14]. 
This means the final Fock space is not equivalent to but completely different from the initial 
Fock space: 

< Oout\ Oin >= . (29) 

This also suggests that the physics completely changes by the emergence of wave-singular naked 
singularities. We can see in the subsequent section that this is indeed the case from the study 
of the vacuum expectation value of the energy-momentum tensor before and after the incident 
at t = 0. 



3.4 Thunderbolts 

Let us see the physical effects of the particle creations by the sudden change of the boundary 
conditions, computing the vacuum expectation value of the energy-momentum tensor. It is 
convenient to introduce the advanced and the retarded time z± :— t ± x, respectively. Then, 
from Eqs. (7) and (8), we can compute the renormalized energy- momentum tensor in t < 
region as 

< in \ T±± | in > t<0 = < in \ (<9±0) 2 | in > 

oo 

n + < in \ {d±hf | in > 

n=l 

h- (30) 



CO 

-Y 



where we have adopted the C-fimction regularization [15]. The right hand side is the Casimir 
energy due to the finiteness of the spatial section. 

On the other hand, from Eq. (23), and the Bogoliubov transformation (24) and (25), we 
can have in t > region, 

1 ( I 1 °° l\ 

<0 m |T ±± |0 m > i>0 = -— + -<0 m |g 2 |0 m >+-V- ){S(z ± )-5(z ± -n)} 2 ,(31) 

\ n=l / 
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with the aid of the formula 



E 

l = — oo 

(i.odd) 



-ilz 



-7TI < Q(z) - Q(Z - 7T) 



E 

/= — oo 
(i.odd) 



-ilz 



71 {S(z) — 5(z — 7r)} 



(32) 



The first term in the right hand side of Eq. (31) is the Casimir energy, which is present also in 
the t < case. From the second term, it turns out that the energy density diverges along the 
null boundaries z + = tt and z_ = 0, which form the Cauchy horizon for the t = initial surface 
as illustrated in Fig. 1. Besides the logarithmically divergent factor, the square of the delta 
function implies that not only the energy density but also the total energy will be unbounded. 
In other words, classical observers passing through the Cauchy horizon will see the thunderbolts. 



t + X =TZ 



t-x = 



t = 




- x 



t<0 



Figure 1: The solid and dashed vertical lines denote the boundaries at 
x = 0, 7T, where the Neumann and the Dirichlet boundary conditions 
are imposed, respectively. The dashed vertical lines mimic the emerging 
naked singularities. Thunderbolts, denoted by zigzag lines, appear along 
the null lines, z + = n and z_ = 0. 



4 Summary and Discussion 

We have studied quantum effects of the sudden change of the boundary conditions which 
mimics the occurrence of wave-singular naked singularities. To be specific, we have dealt 
with the massless scalar field in (1 + l)-dimensional Minkowski spacetime with finite spatial 
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interval. We have shown that the particle creations occur due to the sudden change of the 
boundary condition from the Neumann to the Dirichlet condition at t — 0. Quantum particle 
creation is known to happen when the time-translation operator, or the Hamiltonian is time 
dependent. As we have emphasized, an operator of the spatial part of the wave equation is 
defined as a pair of the action A and its domain of definition T>(A). The sudden change of 
the boundary condition is nothing but the change of the domain T>(A) so that in this sense 
the time-translation operator is time dependent. The behavior of the Bogoliubov coefficients 
in this model implies that the unitarity implementability breaks down, and the Fock space 
for t > is not equivalent to the one for t < 0. Furthermore, having computed the vacuum 
expectation value of the energy-momentum tensor, we explicitly showed that, the singular null 
waves or the thunderbolts appear along the null lines emanating from the edges of the emerging 
singularities, while the energy-momentum tensor is the Casimir energy inside the future domain 
of dependence of the initial hypersurface. 

Our result suggests that the Cauchy horizon associated with wave-singular naked singular- 
ities will be destroyed by the thunderbolts and the global hyperbolicity will be restored by the 
quantum effects. Thus this will be a support for quantum cosmic censorship stating that wave- 
regular naked singularities could appear but wave-singular ones are not admitted to emerge 
when quantum theoretical effects are taken into account. 

Although the present analysis is restricted to the simplest case such that the boundary 
conditions change from the Neumann condition to the Dirichlet condition, it is possible to 
explicitly work out the general boundary condition case. Since U (N) parameters which appear 
in general boundary condition characterize the hair of a wave-singular naked singularity, the 
calculation of the general boundary conditions case will tell us the relation between the strength 
of thunderbolt and a hair of the wave-singular naked singularity. A brief discussion of the general 
boundary condition case is given in the appendix. 

In the present work we have examined only the (1 + l)-dimensional model, for explicit 
demonstration. The results might be applicable even to the 4-dimensional spacetime case, 
following the discussion by Anderson and DeWitt [8]. The point is that the appearance of the 
thunderbolts is due to the sudden change of the domain T>(A) of A at the initial surface. The 
discontinuous change of the basis functions {f n } of T>(A) implies a delta-function like behavior 
of the gradient V/ n at the initial surface. Since the energy-momentum tensor is bilinear in 
V/„, the square of the delta-function, like the second term of Eq. (31), are naturally expected 
to appear in a 4-dimensional setup. 

So far there have appeared many works on quantum effects around a Cauchy horizon, show- 
ing the divergence of the expectation value of the energy-momentum tensor along or near the 
horizon (See e.g., Review [16] and references therein). The particle creation occurs as a result 
of the spacetime dynamics just before the instant of the naked singularity formation by gravi- 
tational collapse (See e.g., Ref. [17]), but is not directly related to the nature of the singularity 
itself. On the other hand, in the present simple analysis, no dynamics of the spacetime metric 
has been considered. The appearance of the thunderbolts may be interpreted as a result of 
the wave scattering off at t = 0+ by the boundary singularities. We can extend our present 
analysis to a gravitational collapse model in (3 + l)-dimensions. Then as commented above, 
we expect a quantum effect, leading to the thunderbolt. This is a new phenomenon, because it 
is not merely a consequence of dynamics which creates naked singularities but directly reflects 
the character of the wave-singular naked singularities. 
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Appendix 

Here we shall consider the general boundary condition case in which the boundary conditions 
at x = 0, 7r are given as bd x <f) = <p before the incident t = but suddenly changes to ad x <f) = 
after t = 0, with a, b being constants. Note that when a, b = 0, the field satisfies the Dirichlet 
boundary condition and when a, b — > oo, the Neumann boundary condition. 

The positive frequency (in)- and (out)- mode functions are respectively given by 

f^ ut \t,x) = —= (sin mx + am cos mx) , in t > , (33) 

^/ixm y/\ + a 2 m 2 

fi in \t, x) = (sinnx + bn cosnx) , in t < , (34) 



so that bd x fn n ^ = fn"^ and ad x jm u€) = /i™'' at x — 0, n. Besides these positive frequency 
modes, when a or b — > oo, there exist zero-modes. But here we do not treat the contributions 
from zero-modes, just for simplicity. 

Then, following the same step in Section 3.2, we can find the Bogoliubov coefficients: 

1 + abn 2 s /o,x 
VI + a l n l \J\ + b l n l 

2 (h \ v 7 ™^ 1 x 
a mn - -{o - a) =-j===== — o m ± n:odd , [6b 

vr VI + a 2 m 2 V 1 + b 2 n 2 m - n 

Pmn - —\P - a) 7= j = = 7= j = = ; ^m±n:odd • (37) 

vr VI + a m 2 V 1 + & 2 n 2 m + n 

From this we can immediately see that, if ab ^ 0, then J^^Li Em=i l/^nm| 2 is finite, while if 
ab = (a 7^ 6), then it is logarithmically divergent, hence in this case the Fock space of the 
(in)-states and that of the (out)-states are unitarily inequivalent. 

The vacuum expectation value of the energy-momentum tensor can be computed in the 
same way in Section 3.4. The result can be written down in terms of Lerch's transcendent, 
which is rather complicated and not so illuminating for drawing the physical consequence. 
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